arXiv:hep-ph/0512329vl 26 Dec 2005 


Measurement of the exclusive light hadron decays of the in e'''e experiments 


P. Wang, X.H. Mo, and C.Z. Yuan 

Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 

(Dated: February 2, 2008) 

The measurement of the exclusive light hadron decays of the if" in e'’'e“ experiments with 
significant interference between the ip" and non-resonance continuum amplitudes is discussed. The 
radiative correction and the Monte Carlo simulation are studied. A possible scheme to verify the 
destructive interference is proposed for the detectors with energy-momentum resolution of (1 ~ 2)%. 

Keywords: ip" exclusive decay; interference; Monte Carlo simulation; e'''e“ experiment. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 


I. INTRODUCTION 


The study of the charmonium has been revived due 
to large data samples collected by CLEOc and BESII as 
well as by the B-factories. A prominent physics which 
has drawn interest for more than two decades is the 
very small branching fractions of pir and other vector- 
pseudoscalar modes in ip' decays compared with their 
large branching fractions in J/ip decays. One proposal to 
solve this puzzle is the 2S-1D states mixing scenario Q 
which predicts enhanced rate of pn mode in ip" decays. 
This scenario is extended to other decay modes besides 
p-K, it foresees possible large partial widthes for the light 
hadron modes in ip" decays [2,12 ■ 

Recently, CLEOc and BESII reported the search for 
various light hadron decays of the ip" ^ |E 12' These 
experiments produce the ip" in e^e“ collision so the ra¬ 
diative correction must be taken into account in the data 
analysis. Another feature of the data analysis of the reso¬ 
nances produced in e^e“ collision is the need to consider 
the non-resonance continuum amplitude. Eor example, 
in the e+e” ^ pir data collected at the ip" mass, if the 
2S-1D mixing scenario gives correct B{ip" pn), then 
the non-resonance continuum amplitude is comparable to 
the Ip" decay amplitude, the measured cross section is the 
result of the interference of the two 0. For the narrow 
resonances like J/ip, ip', T{1S), T{2S), and T(3S'), the 
energy spread of the e“''e“ colliders must also be consid¬ 
ered. But in this work, we concentrate on the wide res¬ 
onances like Ip", of which the width is much wider than 
the hnite energy spread of the e+e“ colliders, so the en¬ 
ergy spread does not change the observed cross section in 
any signihcant way. In this paper, we pay special atten¬ 
tion to the circumstance that the interference between 
the amplitudes of the resonance and continuum has sub¬ 
stantial contribution, particularly the circumstance that 
such interference is destructive. We present the char¬ 
acteristic features of the radiative correction under such 
circumstance. In the forthcoming sections, we begin with 
the amplitudes for the observed processes in e+e" experi¬ 
ments and the parametrization of them. Then we present 
the general properties of radiative correction. Next we 
turn to the Monte Carlo simulation. We discuss the in¬ 
variant mass distribution of the final hadron systems for 


the continuum process and for the pure resonance pro¬ 
cess, as well as for the circumstance in which both the 
resonance and continuum amplitudes exist and there is 
signihcant interference effect between them. Einally, we 
propose a possible scheme to verify the destructive inter¬ 
ference with the data collected at the energy of the ip" 
mass for a detector with the energy-momentum resolu¬ 
tion of (1 ~ 2)%. 

Avoiding complexity but without losing generality, we 
restrict our discussions on two situations in which the 
data are collected either off resonance at continuum, or 
at the energy of the resonance mass. Eor the experimen¬ 
tal setting in which ip" is scanned, the technique details 
are more complicated. We shall leave its study to a future 
work. Through out the paper, the resonance ip" is taken 
as an example, but the analysis can be easily extended to 
other resonances with their width to mass ratios compa¬ 
rable to the energy-momentum resolution of the detector. 


II. THREE AMPLITUDES IN e+e" 
EXPERIMENTS 


The OZI suppressed decays of the ip" into light hadrons 
are via strong and electromagnetic interactions. In gen¬ 
eral, the cross section of e+e” to a certain final state at 
the resonance is expressed in the Born order by 

o-b(s) = ^^^|a3g(s)+ a-y(s)piP(s). (I) 

In the above equation, ,/s is the center-of-mass (C.M.) 
energy of e“''e“, a is the QED fine structure constant, 
a^g{s) and a.y(s) denote the amplitudes in which the reso¬ 
nance decays via strong and electromagnetic interactions 
respectively, V{s) is the phase space factor for the final 
states. However, in e“'"e“ colliding experiments, the con¬ 
tinuum process 

e~'’e~ ^ 7 * ^ hadrons 


may produce the same final hadronic states as the reso¬ 
nance decays do. We denote its amplitude by Uc, then 
the cross section becomes Hi 


o-b(s) 


Attso^ 


|a3g(s) + 07 ( 5 ) + ac(s)P 'P(s). 


( 2 ) 
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FIG. 1: The three classes of diagrams of e'*"e“ ^ 

light hadrons at charmonium resonance. The charmonium 
state is represented by a charm quark loop. 


So what truly contribute to the experimentally measured 
cross section are three classes of diagrams, i.e. the strong 
interaction presumably through three-gluon annihilation, 
the electromagnetic interaction through the annihilation 
of cc pair into a virtual photon, and the one-photon con¬ 
tinuum process, as illustrated in Fig.^ where the charm 
loops stand for the charmonium state. To analyze the ex¬ 
perimental results, we must take these three amplitudes 
into account. 

For an exclusive mode, Oc can be expressed by 

ac{s) = -, (3) 

s 

where J-{s) is the electromagnetic form factor which may 
have a dimension depending on the final states, here we 
define iF{s)V{s) to be dimensionless. We adopt the con¬ 
vention that ttc is real. Since is due to the resonance, 
it is expressed in the Breit-Wigner form 


_ 3FeeJ^(s)/(av^) 


( 4 ) 


where M and F* are the mass and the total width of the 
resonance. Fee is its partial width to e^e“. As a^g is 
also due to the resonance decays, it can be parametrized 
relative to by a complex factor 

C = laag/a^-le*'^, (5) 


as 


039 (s) = C- 


3ree.F(s)/(a7i) 
s- + iMTt 


( 6 ) 


Neglecting double OZI suppressed processes, the am¬ 
plitudes of ijj" decays into pairs of vector-pseudoscalar 
mesons (VP) are parametrized in terms of three pa¬ 
rameters: the strong amplitude g, the electromagnetic 
amplitude e and the SU{3) symmetry breaking factor 
(i-sg) [laiia. These are listed in Table ^ for some 
of the measured decay modes together with the mea¬ 
sured values by CLEOc 0 and BESII II i Simi¬ 
lar parametrization also applies to other OZI suppressed 


two-body decays which conserve the generalized C-parity 
by appropriate change of labelling [^. As for the de¬ 
cays to pairs of pseudoscalar-pseudoscalar mesons {PP) 
which violate the generalized C-parity, the amplitudes 
are parametrized in terms of two parameters: the elec¬ 
tromagnetic arnplitude E and strong SU (3) breaking am¬ 
plitude M [T^ . These are listed in Tab le 1^ to gether 

with the BESII and CLEOc measurements |l2 . ll3| . Sim¬ 
ilar parametrization can be extended to pairs of vector- 
vector mesons (VV). The amplitudes a^g and aj can be 
expressed in terms of this parametrization scheme. For 
example, for pn mode, asg = g, = e, \C\ = |g/e|; for 
K*^KO, asg = g{l-Sg), = -2e, \C\ = |5(l-Sg)/(2e)|; 
while for K'^K~, a^g = = E. \C\ = -^|AI/E|. 

For VP final states, and other final states conserving the 
generalized C-parity, define 

«,=.rg(f), 

while for PP and VV final states, define 



we have (> = 9g if in Tabled the sign between g and e is 
positive; and (j) = 9g + 180° if the sign between g and e is 
negative; while (j) = 9g for K^K~ ■ in Table im To include 
Uc in the formula, simply replace a-y with (a^ + Uc). 

If both the strong and electromagnetic interactions ex¬ 
ist, the partial width of i('' ^ / is calculated by 

Tf = T,,\P{s)\^-\l+C\^V{s). (7) 

If the decay to the final state only goes via electromag¬ 
netic interaction, then 

Tf = T,,\ns)? -r^s). ( 8 ) 

As for the decay to KgKl mode which only goes via 
strong interaction (here we assume that there is no ex¬ 
cited (j) states in the vicinity), iF(s) = 0 since this mode 
does not couple to a virtual photon. Nevertheless, ac¬ 
cording to the parametrization in Table im the decay 
modes K^K~ and have the same strong decay 

amplitude, while tt+tt" and K'^K~ have the same elec¬ 
tromagnetic decay amplitude, which is consistent with 
the recently measured data listed in Table HD within the 
experimental errors Under such assumption, \C\ can 
be determined by the equation 

F^o^o =ree|.F,+ .-(s)MCp7^(s). (9) 

Then with the measurement of K'^K~ mode, (/> is ob¬ 
tained. 

If the data is taken at the energy of the ijj'' mass, i.e. 
s = M^, from Eq. 0) we see that has a phase of 
—90° relative to Uc. So there is no interference between 
ttc and a^. For the decay modes which go only via elec¬ 
tromagnetic interaction, e.g. wtt^, pg, pp' and tt+tt", the 






3 


TABLE I: The experimental results from CLEOc Q and BESII Hi on -i/i" decays into VP modes and biiv. Also listed is 
the parametrization of their amplitudes. Neglecting double OZI suppressed processes, the amplitudes of t/j'' decays into VP 
final states are parametrized by three terms: the strong amplitude g, the electromagnetic amplitude e and the SU(3) breaking 
factor (1 — Sj) |l(l.Il l||. Similar parametrization can be extended to other OZI suppressed two-body decays which conserve the 
generalized C-parity by appropriate change of labelling. The physical states g and g' are expressed by their quark components 
as I 77 ) = Xn^\uu + dd) -b y,,|ss) and \g') = X^i^\uu -b dd) -b y,,'|ss). 


Channel Amplitude CLEOc BESII 





cr(3.671 GeV) [pb] 

cr(3.773 GeV) [pb] 

ct( 3.650 GeV) [pb] 

a(3.773 GeV) [pb] 

VP 

+ 

1 

1 

c 

p tt"'" 

g + e 

8.0+i'4±0.9 

4.4 ±0.3 ±0.5 

< 25 

< 6.0 

UTT® 


3e 

15.2m ± 1.5 

14.6 ± 0.6 ± 1.5 

24.21“ ± 4.3 

10.7m ± 1.7 

CpTV^ 


0 

< 2.2 

< 0.2 



pg 


3eXg 

lO.Oti l ± 1.0 

10.3 ± 0.5 ± 1.0 

8.llimi-l 

7.8l3'5 ± 0.08 

OJT] 


(g + e)Xr, 

2.3t);? ±0.5 

0.4l°;mo.l 



(pg 


[5(1 “ 2sg) — 2e]y,, 

2.1+) ® ± 0.2 

4.5 ±0.5 ±0.5 



pg' 


3eXgi 

2.lttlE0.2 

3.8t^m0-4 

< 89 

< 28 

ujg' 


(g + e)Xgi 

< 17.1 

0.6lomO-6 



(pg' 

OrO 

[5(1 - 2s,) - 2e]y,, 

< 12.6 

2.5t);) ±0.4 



K*°K0, K* 

9{'^-Sg)-2e 

23.5t^m3-l 

23.5 ± 1.1 ±3.1 



K*+K-, K* 

-K+ 

9{^-Sg) + e 

l.Oio.) ± 0.5 

< 0.6 




AP 


g + e 7 . 9 t 2.5 ± 1-8 6.3 ± 0.7 ± 1.5 


TABLE II: The measured ■k~ , K'^K~ cross sections at 
Ecm = 3.671 GeV by CLEOc fl^ and the upper limit of KsKl 
cross sections at Ecm = 3.65 GeV by BES lo . Also listed is 
the parametrization of ip” decays into PP Hnal states. These 
amplitudes are parametrized in terms of electromagnetic am¬ 
plitude E and strong SU{S) breaking amplitude M HD,in 
Similar parametrization can be applied to VV Hnal states. 


final state 

G (3.671 GeV) 

amplitude 

■R^TT~ 

9.0 ± 1.8 ± 1.3 pb 

E 

K+K- 

5.7 ±0.7 ±0.3 pb 

E+^M 

K%Kl 

(< 5.9 pb at 90% G.L. at 3.65 GeV) 



interference between the resonance and continuum can 
be neglected. (Under such circumstance, the interference 
is still non-vanishing due to two reasons: first, for the 
practical reason in the experiments, the data is usually 
taken at where the maximum inclusive hadron cross sec¬ 
tion is, which in gen eral does not coincide with the mass 
of the resonance [ij; second, even if the data is collected 
at the energy of the resonance mass, the interference is 
non-vanishing because of radiative correction. This will 
be proved later on in this paper. For the narrow reso¬ 
nances with widthes smaller than the energy spread of 
the e+e” colliders, the smearing of the C.M. energy also 
results in a non-vanishing interference term. But these 
are beyond the concern for the accuracy of current exper¬ 
iments.) Under such circumstance, in the data analysis 
we simply subtract the continuum cross section from the 
cross section measured on top of the resonance to get the 
resonance cross section. The ratio of the resonance cross 
section of a particular final state to the total resonance 


cross section gives the branching fraction of this mode. 
For the ip”, is very small compared to a^. This is seen 
that if s = M2, 

|a^(M|„)/a,(M2„)| = ^B{iP” ^ e+e”). 

With the measured value of B{ip” e“''e“) = (1.12 ± 
0.17) X10-5 m, 

|a.^(M|„)/ae(M|„)| « 4.6 x IQ-^. 

So can be neglected. For those modes which only 
go via electromagnetic interaction, the measured cross 
sections at the ip” mass almost entirely come from the 
non-resonance continuum amplitude ac- This is demon¬ 
strated by the experimental results on pry, and pg' 
modes in Table 13 where their cross sections measured at 
the Ip” peak are consistent with the ones measured off 
the resonance within experimental errors with the later 
scaled for s dependence. 

But for other final states which have contributions from 
a^g besides a-,., there could be interference between a^g 
and Oc as well as between a^g and a^. Since for the ip”, 
a-y is very small compared to Oc, so only the interfer¬ 
ence between a^g and Oc could be important. Based on 
the analysis of the experimental data, we have suggested 
that the pha se Og is universally —90° in quarkonium de¬ 
cays |Ia[l3- Since at the energy of resonance mass, the 
phase of is —90° relative to Oc, so the relative phase 
between a^g and Oc is either 180° or 0°, depending on 
whether the relative sign between g and e in Tabled or 
between JVl and E in Table HD is plus or minus. The in¬ 
terference between a^g and ac is destructive for the final 
states pn, up, up', K*^K~ + c.c., bi7r, and K~, but 
constructive for <pp, cpp', and K*^K^ + c.c. 
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Destructive interference between asg and Uc means 
that the observed cross section at the energy of the res¬ 
onance mass can be smaller than the continuum cross 
section. The experimental results on pir and WTy modes 
in Tabled demonstrate this interference pattern. 

III. RADIATIVELY CORRECTED CROSS 
SECTION 

The actual description of e“'"e“ to a final hadronic state 
through the annihilation of a virtual photon must in¬ 
corporate radiative correction. Such correction mainly 
comes from the initial state radiation, and for the 
hadronic final state, the final state radiation usually can 


be neglected 0 

The integrated cross section by e“'"e“ collision incorpo¬ 
rating radiative correction is expressed by |l 9|| 

1 —Sm/s 

CTr.c. (s) = J dxF{x, s)aB{s{l — x)) (10) 

0 

where a 3 ( 3 ) is the Born order cross section by Eq. @, 
and in the upper limit of the integration ,ysm is the cut¬ 
off invariant mass of the final state hadron system after 
losing energy to photon emission. F{x, s) is calculated 
to an accuracy of 0.1% in Ref. [T^ : 


F{x,s) = f3x^ ^ 

1 


1 + iP ”1-( -n 


^ ' 32 12 


%i-f) 


+ 0/3 


4(2 — a;) In- 

X 


1 1-1- 3(1 — xY 


ln(l — x) — 6 -I- a; 


( 11 ) 


with 



TT 



( 12 ) 


In the above equations, me is the mass of the electron. 
Here the expression F{x,s) includes the bremsstrahlung 
of an e“'"e“ pair from the initial e“'"e“ state. For x ^ 0, 
Xy/sl2 is approximately equal to the energy carried away 
by the radiated photons. But for x ~ 1, this meaning is 


not valid in the order. The effects of vacuum po¬ 
larization are not included explicitly in Eq. m- Here 
we follow the convention that for hadronic final state, 
the vacuum polarization by leptons and hadrons, includ¬ 
ing vector-meson resonances, is taken into account in the 
form factor F{s) |^ . 

In Eq. (Cni), F{x,s) is positive definite. If s = M^, 
for pure electromagnetic processes, the interference term 
between a-y and Uc in (Tb(s(1 — x)), i.e. 


23? [a.y(s(l — x))ac(s(l — x))] 


23? 

J^(s(l-x)) 

s(l — x) 

2 [s 

(1 - x) - 


3reeJ^(s(l - x))/{a^s{l - x)) 
s{l - x) - + iMTt 

_ 3|.F(S(1-X))|"ree _ 

a[s(l — x)]^/^ [(s(l — x) — + M^Tj] ’ 


(13) 


vanishes only at x = 0, but it is negative elsewhere in the 
integration interval from 0 to 1 — Sm/s. So even in the 
pure electromagnetic processes, with radiative correction, 
the interference term gets a negative value. 

In Eqs. ®,0 and ®, the form factor F{s) is adopted 
to describe the hadronic interaction. Both the Monte 
Carlo simulation and the calculation of the radiatively 
corrected cross section require the knowledge of the form 
factor in the energy range from ^/s down to -y/FY- In 
principle, y/FY can be as low as the production threshold. 


But this requires the input values of the form factor from 
y/s to the production threshold. In Eq. (trnil . F(x, s) has 
been calculated to an accuracy of 0.1% by QED, but for 
virtually all the hadronic final states, we have neither 
the theoretical models nor sufficient experimental data 
to describe the form factors to such high precision. To 
reduce the systematic uncertainty from the Monte Carlo 
simulation, one strategy is to take the value of Sm as close 
to s as possible for the actual event selection criteria, i.e. 
we generate the event sample with Sm just lower than 
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that which are to be selected. In this way, we need only 
to describe the form factor precisely from ^/s down to a 
much higher instead of to the production threshold. 

The variation of the form factor in a small energy in¬ 
terval usually can well be approximated by 

( 14 ) 

with k either derived from theoretical models, or obtained 
from fitting the experimental data at nearby energy. Of 
course one may use a more complicated function other 
than Eq. o for the approximation. 

IV. MONTE CARLO SIMULATION 
A. Soft and hard photon events 

In realization of Monte Carlo simulation, an auxiliary 
parameter xq is introduced to separate the integration 
interval of Eq. m into two parts, (0,a;o) and (ccq,! — 
Sm/s), with xo a small but nonzero value: 

XQ 

O'so/t(s) = J dxF{x,s)aB{s{l - x)) (15) 

0 

and 

1 /<5 

Uhard{s) = J dxF{x,s)aB{s{l - x)) (16) 

Xo 

In the above two equations, all functions and variables 
have the same meaning as in Eq. cni. We have 

The two terms in the above expression are usually called 
soft photon and hard photon cross sections respectively. 
The Monte Carlo program generates the soft photon 
events and hard photon events according to their pro¬ 
portions in the total cross section. In a soft photon 
event, the four-momentum of the radiated photons are 
neglected, and the photons are not generated. Experi¬ 
mentally, this means that in the soft photon events, the 
photons can not be detected. Since in Eq. OS , xq is 
small, so XQ^fsj^l can be identified as the maximum en¬ 
ergy carried away by the undetected photons. This re¬ 
quires that XQyfsj^ must be smaller than the minimum 
energy of the photon which can be detected in the exper¬ 
iment. In a soft photon event, the energy-momentum is 
conserved between the incoming e“'"e“ pair and the final 
state hadron system, so xq must also be smaller than the 
energy-momentum resolution of the detector. Usually xq 
is assigned a nonzero value smaller than 0.01. The out¬ 
come of the Monte Carlo simulation does not depend on 
Xq . In a hard photon event, the four-momentum of the ra¬ 
diated photons are generated, and the energy-momentum 
is conserved with the inclusion of these photons. 


The differential cross sections with the emission of pho¬ 
tons for exclusive hadronic processes are calculated in 
Refs, the calculation of Ref. can also be 

used to generate the hard photon events while the gen¬ 
eration of soft photon events is identical to the genera¬ 
tion of events at Born order. (The inclusive process is 
treated somewhat differently in Ref. [l^.) The authors 
of these references also provide Monte Carlo programs 
based on their calculations. The program based on the 
calculation in Ref. [^, BABAYAG A, g enerates tt+tt" 
events; the program based on Ref. jU, MCGPJ with 
the current version, generates tt+tt”, K~^K~, and KsKb 
events; while the program based on Ref. generates 
hard photon events for tt+tt”, Tr+Tr^Tr*^, 7r+7r“7r+7r“, and 
7 r+7r“7r°7r° processes. These programs achieve the preci¬ 
sion of (0.1 ~ 0.2)%. 

In order to generate events with the presence of both 
the resonance and continuum, some replacements are 
to be made in the fore-mentioned programs. The Born 
order cross section by Eq. must be substituted to 
generate the correct distribution of the invariant mass 
of the final state hadron systems. Although the origi¬ 
nal programs only generate a few hadronic final states, 
the programs based on Refs. mm are written in the 
form that more final states can be added in a straightfor¬ 
ward way. To do this, one need to put the corresponding 
hadronic tensor into the program, where 

= n^x nt, (17) 

with the current of virtual photon transition to the 
final hadron state /. For reference we include its forms 
for some final hadronic states in the appendix. 

B. Distribution of invariant mass of hadrons 

The distribution of the invariant mass of the final 
hadron system Miny is very different between resonance 
and continuum. More profound distinctive feature is in 
the circumstance that there is interference, particularly 
destructive one, between the resonance and continuum. 
In this section, we discuss this distribution. 

1. Resonance 

If the final state does not couple to a virtual photon, 
e.g. KgK'l, then in Eq. |(2J), there is only the a^g term, 
the Born order cross section is expressed by the Breit- 
Wigner formula: 

, . _ 127rreer/ 

'^’'“('®)“ (s-M2)2-HM2r2’ ^ ^ 

where T/ is the partial width to the final state /. 

For narrow resonances, by narrow we mean Tt/M <C 
Xo, Eq. m behaves almost like a 6 function. If the data 
is taken at the energy of the resonance mass, i.e. s = 
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in Eq. Cl, crres(s(l — x)) is very large at a: = 0 and 
very small elsewhere. So in Eq. d, the contribution 
of (Tb(s(1 — x) to the radiatively corrected cross section 
comes solely from a small interval around a; = 0. Sub¬ 
stitute ares(s(l — x)) for aB{s{l — x)) in Eqs. and 
d, asoft is nonzero while ahard virtually vanishes due 
to their integration intervals respectively. This means 
that in the Monte Carlo simulation of a narrow reso¬ 
nance, only the soft photon events are generated. The 
sole effect of the radiative correction is the reduction of 
its height. This will be discussed later in this paper. 

The resonance ip” is not very narrow in this sense. If 
we take xq — 0.01 and Sm = 0.8M^, soft photon events 
are 98.6% of the total. 


2. Non-resonance continuum 

For continuum process, there is a considerable hard 
photon cross section. Take an example of EP final states. 


with xq = 0.01 and Sm = 0.8s, and assuming that the 
form factor varies as a function of energy according to 1/s 
(fc = 1 in Eq. C3)i in the simulation the hard photon 
events are 22% of the total. 


3. Interference between resonance and continuum 


If both the resonance and continuum amplitudes exist, 
the Born order cross section is obtained from Eq. © with 
ttc, a-f and a^g by Eqs. @, 0, and ©■ 


J 


ctb(s(1-x)) = 


47 rs(l — x)a^ 


|ac(s(l - a;)) -I- a 3 g(s(l - x)) -f Ot,(s( 1 - a;))| P(s(l - x)) 


= 2 '^'S (1 — x)a^\T{s{l — a;))pP(s(l — x)) x 


s(l — x) 


+ (C + 1) 


ST, 


e/ ^a-\/s(l - a:)j 


s{l - x) - + iMTt 


■ ( 19 ) 


We pay special attention to the destructive interference 
and take as an extreme exmaple the circumstance that 
a^g and Uc almost cancel to each other in the Born order. 
This happens if 

</> = -90°. (20) 

ee 

Under such circumstance, the partial width of this final 
state is expressed by 

where acon{M‘^), scaled for s dependence to s = M^, is 
the Born order cross section for the non-resonance con¬ 
tinuum. 

To calculate the proportion of soft and hard photon 
events, Eq. is to be substituted into Eqs. 03 and 
dm respectively. Notice that these two equations only 
differ by their integration intervals. For the discussions 
on the Ip" in this work, we take xq = 0.01 which is greater 
than Tt/M. If the resonance amplitude satisfies Eq. (12011 . 
then for Eq. 113, in the integration interval (0 , xq), there 
is almost complete cancellation between the two terms of 
tTB(s(l — a;)); while for Eq. II16II . in the integration inter¬ 
val {xo, 1 — Sm/s), the magnitude of the resonance which 


is the second term of (Tb(s(1 — x)), virtually vanishes. 
This means that the resonance and its interference with 
continuum affect mainly the soft photon cross section; 
while the hard photon cross section is predominately due 
to Uc- In this way, the interference changes the propor¬ 
tions of the soft and hard photon events as well as the 
total radiatively corrected cross section. The destruc¬ 
tive interference between resonance and non-resonance 
continuum reduces the soft photon cross section, which 
means a smaller proportion of soft photon events. 

To illustrate the above discussion quantitatively, we 
take the example that the resonance amplitude satishes 
Eq. 123. Under such circumstance, the proportion of the 
soft photon events is 15% for xg = 0.01 and Sm = 0.8s. 
(This is not the lowest possible proportion of the soft 
photon events. With (p = —90°, the complete cancella¬ 
tion happens only between a^g and Oc. There is still a 
small but non-vanishing a-y left. If a^g cancels out not 
only Oc, but also a.y, the proportion of soft photon events 
can be even smaller.) This means that for a detector 
with energy-momentum resolution xg/2 = 0.5%, among 
the data taken at non-resonance continuum, 78% of the 
events have invariant mass equal to the C.M. energy of 
the incoming while among the data taken at the 

Ip" mass, we may find that the majority of the events 
(85%) have invariant mass smaller than the C.M. en- 
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ergy of the incoming e+e“. If this happens, it indicates 
destructive interference between the resonance and non¬ 
resonance continuum amplitudes. 

Fig.[21shows the probability density as a function of the 
squared invariant mass of the final hadron system . 
For continuum, as well as for no interference or construc¬ 
tive interference between the resonance and continuum 
cases, the maximum probability density occurs (actually 
diverges) at = s, which corresponds to ^ 0 with 
the energy of the emitted photons; but if the destruc¬ 
tive interference between a^g and Oc exists, and it leads 
to almost complete cancellation of the two amplitudes, 
the probability density may have a minimum point near 



FIG. 2: The distribution of the probability density as a func¬ 
tion of Mf^^/s with Minv the invariant mass for a. VP final 
state. The solid line is for continuum. The dashed line is an 
example of destructive interference between asg and with 
the two amplitudes satisfying Eq. iCT . The probability is 
normalized to with Sm ~ 0.8s. 


In the circumstance that at s = the three terms 
in Eq. completely cancel, i.e. a^g -I- a.y -I- Oc = 0, so 
the Born order cross section vanishes, the radiatively cor¬ 
rected cross section ar.c.{s) still gets a nonzero value, as 
long as Sm < s. This is because that in Eq. (tTUll . F[x,s) 
is a positive definite function, while crB(s(l —a;)) vanishes 
only at a; = 0, but remains positive definite elsewhere 
in the integration interval. This leads to the phenom¬ 
ena in e“'"e“ experiments: a detector with finite energy- 
momentum resolution always observes a non-vanishing 
cross section, even if the Born order cross section van¬ 
ishes due to destructive interference between the reso¬ 
nance and continuum. 


C. The calculated efficiency 

According to Eq. dia, CTr.c. is a function of Sm as well 
as s. For the reason which we discussed in Sec. im in the 
Monte Carlo simulation, usually we generate events to a 


cut-off invariant mass y/Sm, so the calculated efficiency e 
is also a function of Sm- 

If the observed number of event is denoted as N, cor¬ 
responding integrated luminosity is denoted as £, then 
we have 


^m) ‘ 


here and in the following discussions, we explicitly in¬ 
dicate the dependence of Ur.c. and e on Sm- The above 
equation can also be expressed as 

N 

— ^r.c.(s,577i) * 6(Sm) ■ 

The left side of the above expression is an experimentally 
measured quantity. The product crr,c.(Sj ^m) ‘ e(sm) does 
not depend on Sm- 

Very often, the experimental results are presented in 
terms of Born order cross section, particularly the off 
resonance continuum cross section. For the data taken 
off the resonance, we measure the electromagnetic form 
factor of the final state, which is simply related to the 
Born order cross section. A so-called radiative correction 
factor is introduced as 

f f \ <^r.cXs,Sm) /'oo^ 

fiSR{s,Sm) = -rr—• (23) 

o-b(s) 

Here if the form factor takes the form of Eq. mi), then 
fisnis, Sm) can be calculated. It does not depend on 
lF(s), although it still depends on k in Eq. itT^ . With 
Eqs. (1221) and (P|l. we have 


^ Born{^^ 


^r.c. ) 

fisnis, Sm) 


N 

£ • e(sm) • fisnis j ^m) 


In the denominator, the product of e(sm) ■ fisnis, Sm) 
cancels out the dependence on Sm- 

The interference between the resonance and continuum 
amplitudes may change the factor fisR in a profound 
way. For non-resonance continuum, with ^/s well above 
the production threshold, if the form factor goes down 
rapidly as s increases (e.g. fc > 1 in Eq. m) , as Sm 
approaches the threshold, f/sR is usually greater than 1; 
but if Sm is taken close to s, then fjsR can always be 
smaller than 1, since as Sm —>■ s, crr.c.(s, Sm) —>■ 0. Eor 
a resonance, frs R do es not depend on k. It is roughly 
approximated by [3 


fiSR 



. 3 a 

1 -b -/3 H— 
4 TT 




Here (3 is given by Eq. mi- If Ft <C M, the value of 
this expression is less than 1. It means that the ini¬ 
tial state radiation reduces the height of the resonance. 
If there is significant interference between the resonance 
and continuum, fisR may take any value between the 
order of 1 to infinity. It becomes infinity in the circum¬ 
stance that (7 b (s) vanishes due to destructive interference 


( 24 ) 










between the resonance and continuum. As discussed in 
Sec. II V o'r.c.(s, Sm) still gets a nonzero value in this 
circumstance. 

For illustrative purpose, here we take an example that 
in the Monte Carlo simulation of e+e” collision at ^/s = 
M^ii, the VP final state events with the cut-off invariant 
mass of 0.9M^ii are generated, and the form factor varies 
as a function of the energy according to 1/s, while in 
the data selection, the invariant mass of the hnal VP 
particles is required to be greater than 98%^/s. Under 
these conditions, for the continuum cross section, 87.4% 
of the generated events survives and fisR = 0.946; for 
the 4’" resonance, 99.9% of the generated events is left 
and fisR = 0.716; while for the destructive interference 
between a^g and Uc which satisfies Eq. (jSnj, only 49.0% 
of the generated events have invariant mass greater than 
98%My,», but fisR = 65.5. 


V. MEASUREMENT IN THE PRESENCE OF 
INTERFERENCE 

The above discussions lead to a profound feature of 
the experimental measurement in the presence of inter¬ 
ference between the resonance and non-resonance contin¬ 
uum: the complete determination of the branching frac¬ 
tion must come together with the determination of the 
phase between the resonance and continuum by scanned 
data around the resonance peak. The data must be 
taken at least at four energy points, because there are 
three quantities which must be determined simultane¬ 
ously: \C\ and </. At the same time, the form of 

dependence of the observed cross section on \C\ is quadric. 
Herein the non-resonance continuum amplitude or equiv¬ 
alently iF(M^) is determined by the data taken at con¬ 
tinuum. In the treatment of the data taken around the 
resonance, if both strong and electromagnetic interac¬ 
tions exist, the Monte Carlo generator requires the input 
of \C\ and 0, so the data analysis is an iterative process. 
The usual procedure is to fix (/, and varying |C|, until 
the calculated efficiency e(sm) by the Monte Carlo and 
the radiatively corrected cross section ctj-.c. (s, Sm) satisfy 
Eq. 122). Then the partial width is obtained by Eq. 0. 

If the data are taken at two points, i.e. one at contin¬ 
uum off the resonance and the other one at the energy of 
the resonance mass, only a relation between \C \ and (j) can 
be obtained. The solution of \C\ and 4> is differentiated 
into two circumstances, depending on the relative mag¬ 
nitudes of the observed total cross section at and contin¬ 
uum cross section a^. Here CTc is the radiatively corrected 
cross section of the non-resonance continuum calculated 
with the form factor at the energy ^/s = M^". As in 
the Born order cross sections, if at > ac, 4> can take any 
value from —180° to 180°, and for each value of 4 there 
is one and only one solution for |C|; on the other hand, if 
at < CTci 4 is constrained within a range around —90°, in 
which every possible value of 4 corresponds to two solu¬ 
tions of |C|. A formal discussion is left into the appendix. 


In such measurement, \C\ is determined versus the phase 

in a two dimensional curve. The recommended way is 
to start from 4 = —90°, since there is always at least one 
solution of \C\. If the obtained at is smaller than CTcj the 
second solution must be searched. Then for the 4 val¬ 
ues greater and smaller than —90°, the solntions of \C\ 
are found similarly. Thus the curve relating \C\ with 4 is 
obtained point by point. 

Although (Tr.c.(s,Sm) and e{sm) depend on \C\ and </, 
for different solutions of \C\ and 4 which fit the exper¬ 
imental data, ar.c.{s, Sm) and e{sm) only depend on \C\ 
and 4 weakly. This can be understood by noticing that 
the efficiency can be expressed in terms of the efficiency 
for soft photon events esoftixo), and the one for hard 
photon events Chardism)- Apparently we have 

^r.c. (s, — ^softi^Td^o)^soft{^o)P^hard{,^ 

n—rr-. 

As discussed in Sec. if we take ccp > Ft /M, the dis¬ 

tribution of the hard photon events and the hard photon 
cross section are predominately due to Oc, so ahard{s, Sm) 
and Chardism) are almost independent of \C\ and 4 - As for 
the soft photon events, the distribution follows the Born 
order differential cross section, so eso/t(a^o) does not de¬ 
pend on \C\ and 4 either. So in Eq. 12511 . for a rough 
approximation, only asoftis,xo) depends on \C\ and 4 - 
The solutions of |C| and 4 are to satisfy the equation 

N 

— ^soft{s)^soft P ^hardi^: S7n)^hard{^m) ■ 

Here N/C is given by the experiment, only asoft{s,xo) 
depends on \C\ and 0, so asoft{s,xo) remains as a con¬ 
stant for any possible solution of \C\ and 4 which fit the 
data. Since (^5 ^m) — ^soft(^^j^o)P^hard(^^^^m)i and 
N/C = ar.c.{s,Sm)eiSm), cTr.cXs.Sm) and e{Sm) are al¬ 
most constants as well. This property helps to find other 
solutions with different \C\ and 4 values once one of the 
solution is found. It also makes the iterative process con¬ 
verge very fast. 

VI. VERIFY THE INTERFERENCE WITH THE 
DATA AT 4” 

The discussions in Sec. EYH leads to a scheme which 
could verify the destructive interference between 4 ” and 
non-resonance continuum with only the data at 4 " peak. 
Modern detectors with a CsI(Tl) calorimeter and a mag¬ 
netic field of 1 Tesla or more, like CLEOc and BESHI, 
measure the energy-momentum with resolution of 1%, 
which is comparable to the ratio F^/z/M^//. In this 
scheme the invariant mass distribution of the hadrons 
is measured, and in the event selection, requirement of 
the invariant mass, Minv, greater than a certain value, 
Mcut, is applied. For non-resonance continuum, as Mcut 
is loosed from 0.99'\/s to, e.g. 0.95-\/s or 0.90-^5, the 
number of events increases slowly; for no interference or 
constructive interference between the resonance and con¬ 
tinuum, the number of events increases even slower; for 
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TABLE III: The variation of the cross section as a function 
of Mcut for the non-resonance continuum and the destructive 
interference between azg and Oc with the two amplitudes sat¬ 
isfying Eq. I121)L taking the cross sections with Mcut = 0.99-y% 
as the unit. 

Mcutl-Js 

destructive interference 

continuum 

0.99 

1.00 


1.00 

0.98 

1.61 


1.06 

0.97 

2.00 


1.09 

0.96 

2.30 


1.12 

0.95 

2.53 


1.14 

0.90 

3.33 


1.21 


pure resonance, the number of events does not increase 
at all; on the other hand, if the destructive interference 
leads to substantial cancellation between the resonance 
and continuum amplitudes, as Mcut is lowered, the num¬ 
ber of events increases rapidly. Table nm gives the the 
cross section as a function of Mcut, taking the cross sec¬ 
tion with Mcut = 0.99-v/i as the unit. Listed cross sec¬ 
tions are due to the non-resonance continuum and the 
destructive interference between 03 ^ and Oc with the am¬ 
plitudes satisfying Eq. CT . It is assumed that the final 
state is y P and the form factor varies as a function of en¬ 
ergy according to 1/s. Fig.|21shows the cross section as a 
function of Mcut for these two circumstances, normalized 
to the cross section with Mcut = 0.99-v/s. In Fig. Eland 
Table IHTl we see that, as Mcut is loosed from 99%y^ to 
90%-\/s, the cross section of the continuum increases by 
merely 21 %; but for the destructive interference, it could 
increase by more than 3 times. In the experimental data, 
if the number of events increases rapidly as one lowers 
Mcut in event selection, it indicates destructive interfer¬ 
ence between tit" and continuum with the two amplitudes 
in comparable strength. 


.0 3.5 

0 ) 3 

<0 

58 2.5 

o ^ 

^ 1.5 


TO 

g 

o 

c 


1 - 
0.5 f 


0?99 


0.94 

cut c.m. 


0.90 


FIG. 3: The cross section as a function of Mcut normalized to 
the cross section with Mcut = 0.99y%. Solid line is for contin¬ 
uum. Dashed line is for destructive interference between a^g 
and Oc with the two amplitudes satisfying Eq. 11201 . 


To simulate the experimental situation, in Fig. El the 
interval of the invariant mass of the final state hadron 
system from 100% to 90% of the e+e” C.M. energy is 
divided into 20 bins, and the probabilities of the hadrons 
in the bins are plotted for the events at continuum and 
at the energy of ip" mass with the destructive interfer¬ 
ence between a^g and Oc satisfying Eq. it^ . Here off 
resonance at continuum, the events are highly concen¬ 
trated in the last bin with Minv/Ec.m. = 100 %; while for 
destructive interference, the events are distributed more 
flatly among the bins. If there is no interference or if the 
interference is constructive, the events are even more con¬ 
centrated in the last bin than off resonance at continuum 
in the plot. 
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0.9 0.92 0.94 0.96 0.98 

M „ 

inv c.m. 


FIG. 4: The interval of the invariant mass of the final state 
hadron system from 100% to 90% of the e''"e“ C.M. energy 
is divided into 20 bins, the probabilities of the hadron event 
in each bin are plotted for continuum (dots) and for the tp” 
with the destructive interference between a^g and ac satisfying 
Eq. II 2 OII (boxes). 


This scheme requires the selected data sample be free 
from background contamination. The most important 
background comes from the radiative tail of the tp'. For 
the data taken at the tp" peak, the radiative tail due to 
the tp' has an invariant mass which is 97.8% of the e+e” 
C.M. energy. The radiative tail of the tp' at the energy 
of the tp" mass has a total cross section of 2.6 nb, so this 
scheme is feasible for those decay modes like pir, K*~^K~, 
and utr] with branching fractions in tp' decays no more 
than the order of 10 “® since the background of 

these modes from radiative tail of the tp' is at the order of 
0.1 pb. But for those modes which have large branching 
fractions in tp' decays, e.g. 6 i 7 r, the contribution of the 
tp' tail must be considered carefully. 


VII. SUMMERY 

In this paper, we examined the radiative correction and 
Monte Carlo simulation for the measurement of tp" exclu¬ 
sive decays into light hadrons in e“'"e“ experiments. We 
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draw special attention on the interference effect between 
the tp" resonance and non-resonance continuum ampli¬ 
tudes. We analyzed how the interference, particularly the 
destructive interference may change the invariant mass 
distribution of the final state hadrons. The discussions 
also lead to a possible scheme to verify the destructive 
interference using only the data taken at the energy of 
the Ip" mass for the detectors with the energy-momentum 
resolution of (1 ~ 2)%. We suggest this scheme be ap¬ 
plied on decays such as lot], ujrj', K*~^K~, and /ott by 
future BESIII experiment. 

APPENDIX A: THE HADRONIC CURRENT 

The current of virtual photon transition to the final 
state / is defined as 

= {f\outJr\^). 

For reference we include its forms for some final hadronic 
states in the following. A complete list of for all 
possible two-body final states can be found in Refs. 


1. Pseudoscalar-pseudoscalar 

For the final states tt+tt” or K^K~, 

T-if, = J^p{s){p+ 


where p+ and p- are the four-momentum vectors of tt’’' 
(Ar+) and 7r“ {K~) respectively, and ^p{s) is the tt or 
K form factor at the energy scale s = {p+ +p-)‘^. 


2. Vector-pseudoscalar 

For the VP final states, 

“Hfi — ^Vp{,^')^pal3^PYPp^ ^ 

with the completely antisymmetric unit tensor of 

fourth rank, pv and pp are the four-momentum vectors 
of the vector and pseudoscalar mesons respectively, e is 
the polarization of the vector meson, !Fvp{s) is the form 
factor at the energy scale s = {pv +pp)^. 

3. 7r''"7r~7r'’ with pir intermediate states 

For the 3-body final states 7r+7r“7r°, 

~ ^37r{^')^paf3^P-\-P—PQ: 

with Pp, p_ and po the four-momentum vectors of 7r+, 
7r“ and respectively, T^Tr{s) the form factor at the 
energy scale s = (p+ + p- + po)^. 

We may also include the three intermediate states 
p'^n~, p~TT'^, p^TT^ and their interference into by mul¬ 
tiplying a factor 


-b 


m„ 


+ 


m„ 


pX^ - ml + iTp{p\_)s/mp ' p% - -b irp(p%)s/mp ' p'io - mf, + iTp(p'ig)s/r 


where irip and rp(s) are the mass and energy-dependent 
width of p, p+- = P++P-, p+o = p+ -b Po, and p_o = 
P- +Po- 


APPENDIX B: SOLUTIONS OF \C\ AND 0 

In this section, we discuss the possible solutions of \C\ 
and (p and so the branching fraction of the resonance 
decays if the experimental data are available at only two 
energies, one off the resonance and the other one at the 
energy of resonance mass. We begin with Fq. and 
define 


Here it is more convenient to consider instead of \C\ and 
(p, the total resonance amplitude by 

p = |C -b 1|; 

Ip = arg(C -b 1) . 

Here p = |(a 3 p -b a.y)/a.y| is the total resonance ampli¬ 
tude normalized to the electromagnetic amplitude, and 
Ip = arg((a 3 g-ba-),)/a-),) is the phase of the total resonance 
amplitude relative to the electromagnetic amplitude. As 
discussed in Sec. m for ip", ja-yl <C locI, so can be ne¬ 
glected, only the amplitudes of a^g and Qc are important. 
These two amplitude have significant interference if their 
magnitude are comparable, i.e. 


U = 


s(l — x) 


+ (C + 1) 


3r, 




s(l - x) - M2 -b iMTt 


(HI) 


1^1 

For the discussion on the interference, we only need to 
consider the circumstance that \C\ 3> 1. Then we have 
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approximately p k, \C\ and '0 « 0. This means that 
since |a^| jad, if |a 3 g| is comparable with jad, then 
osg + a-y ~ osg, i.c. the total resonance amplitude is 
approximately equal to the amplitude via strong decay. 

For briefness, we introduce the following notations: 


and 


tan A = 


B 

A ■ 


Introduce a variable ^ and let 


1 , 3r 

<7 = -r, -* = 


s(l — x) s(l — x) (B2) 

a = s(l - cc) - b = MTt, 


then we have 


e = q+(\C\e^^+ 1). 


q + ■ 


a + ib 

2 


a + ib 

Eq. (IB3II can be rewritten as 

+ 2qRcos{ip — X) + 

with 

R = - 

tan A = 


(B3) 


(B4) 


Va2 + 62 
b 


To take the radiative correction, we introduce an inte¬ 
gral operator defined as 

1 — Sm/s 

jiG^ j 4„(l-»)a=|.R(.(l-x))PP(.(l-x)). 


then the radiatively corrected cross section becomes 

^ (s) = J dGt^ = . 


(B5) 


Notice 


qRcos{tp — X) < qR = 


qk 


Va2 + 62 ’ 


and use the relation 


dGqR 


< / dGq^ ■ / dG 


+ b‘^ ’ 


we have 


\A42+^ cos(V’ - A) < QK, (B6) 


where 


aqk 


^ J ^ ^ j 


bqk 




= / dG- 


= / dGq^ , 


- cos(d-A), (B7) 

QK 

we obtain an expression similar to Eq. 1B4II : 

+ 2pQK cos i + p^K'^ . (B8) 

In virtue of Eq. iITTtIi . the angle ^ does not have an appar¬ 
ent physical meaning but has a rather complex relation 
to the angle ip. 




FIG. 5: Possible solutions of p and ip for T > Q (a.) and for 
T < Q (b). Here T, Q, pK denote the total, continuum and 
resonance amplitudes, respectively. In the figures ^' = ip' — A, 
where ip' = —ip according to the angle definition between two 
vetors. In (b) the dashed line corresponds to one solution case 
for special . 

Introduce an angle C defined as ^ = tt— then Eq. iTBSl) 
becomes. 


= Q'^ - 2pQK cos C + p^K^ , (B9) 

in which the three variables T, Q and pK form a tri¬ 
angle and obey the law of cosines. If T > Q, C 
take any value between —180° ^ 180°, there is one and 
only one solution for a given value of while if T < Q, 
there are two solutions for a value of Cj but the range 
of C is constrained by sinC < T/Q. Eor sinC = Q/T, 
the two solutions coincide. In Figs. EJa) and we 
show schematically the two cases for T > Q and T < Q 
respectively. 

Next we consider the variation of €{sm)- As discussed 
in Sec. El to determine the efficiency by the Monte Carlo 
simulation, the resonance parameter p and ip are input 
parameters which are also the quantities to be measured. 
So the data analysis is an iterative process. Here we are 
to prove that if the efficiency remains stable for different 
solutions of p and ip, or with only weak dependence on 
these two variables, the above discussion is still valid. 

In Eq. (IB911 . is the total observed cross section, 

which is obtained by experimentally measured quantities 
as 



-b 


C-e 
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where the efficiency e is a function of p and Tp. With¬ 
out losing generality, we only write its dependence on p 
explicitly, i.e. e = e{p). As shown in Sec. 0 e{p) re¬ 
mains stable for different solutions of p and tp, and the 
dependence on these variables is weak, so we take Taylor 
expansion of e in terms of p, and neglect the higher order 
terms: 


We get 


P = 


{Q cos (- 6at)± \J (QcosC - PutY + 


with 


Q2)J jK, 
(B12) 


e(p) 


e(Po) + 


dp 


i.P - Po) + 0[{p - pof] 

P=Po 


^o + VP = eo 




Consider the quantity under the radical 


(B13) 


where 


eo = e(po) - 


de{p) 


dp 


Po 


P=Po 


V = 


de{p) 


dp 


P=Po 


Here the correction term to eo is small, i.e. 


^«1. 

eo 


z = {Q cos C - 6at)'^ + {at - Q'^) 

= 6‘^a^ — {2SQ cos ( — l)at — sin^ ( . 


In order for Eq. (IHT 2 I 1 to have solutions, it requires z > 0. 
For at > this is always true. Notice that p must be 
greater than 0 by definition, in this circumstance, p has 
one and only one solution for any given value of C, with 
no constraint on (p. If at < Q^, in order to have z > 0, 

(BIO) 

we need the condition 


Therefore we have 

N 


£ • eo ■ ( I -l- 


VP 

eo 


N 

C ■ eo 



= at-l 1 


where we define 


CTt 


N 

’ 



CTt > 

5 


{2SQ cos C — I) -|- {26Q cos C — 1)^ + 4(5^Q^ sin^ ( 

2^2 


which imposes constraint on ({. In such case, there are 
two solutions for each allowed value of (p. In Eq. fim . 
rj/eo is small, so is 6 . In the limit <5^0, we find p has 
two solutions when sin^ (p < at < ■ If z = 0, which 

means 


with at the total cross section calculated by the efficiency 
eg. is the continuum cross section which does not 
depend on p and ip. Substituting the above expression of 
back into Eq. (IB9II we have 


CTt 


{2SQ cos C — 1) -|- Y^^(25Q'cos'c^^T)M-4^2Q2'^i^^^ 
2^2 


eo 


K^p^- 2QKcosC-at- ] p + {Q^ - at) = 0 . (Bll) 


the two solutions of Eq. (TmUl becomes one. 
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